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ABSTRACT 

The central region of the galaxy Henize 2-10 hosts a black hole (BH) candidate with a mass 
Log = 6.3 ± 1.1. While this putative black hole does not appear to coincide with any 

central stellar overdensity, it is surrounded by 11 young massive clusters with masses above 10® Mq. 

The availability of high quality data on the structure of the galaxy and the age and mass of the clus¬ 
ters provides excellent initial conditions for studying the dynamical evolution of Henize 2-10’s nucleus. 

Here we present a set of 7V-body simulations in which we model the future evolution of the central 
clusters and the black hole to understand whether and how they will merge to form a nuclear star 
cluster. Nuclear star clusters (NSCs) are present in a majority of galaxies with stellar mass similar 
to Henize 2-10. While the results depend on the choice of initial conditions, we find that a NSC with 
mass Mnsc — 4 — 6 x 10® Mq and effective radius tnsc — 2.6 — 4.1 pc will form within 0.2 Gyr. 

This work is the first showing, in a realistic realization of the host galaxy and its star cluster system, 
that the formation of a bright nucleus is a process that can happen after the formation of a central 
massive BH leading to a composite NSC-I-BH central system. The cluster merging process does not 
significantly affect the kinematics of the BH; when a stationary state is reached its position changes 
by < 1 pc and its velocity by < 2 km s“*. 

Subject headings: galaxies: individual (Henize 2-10), galaxies: nuclear star clusters, galaxies: star 
clusters; methods: numerical. 

1. INTRODUCTION 

The majority of the observed galaxies with luminosi¬ 
ties up to 10**Lq host bright stellar nuclei usually re¬ 
ferred to as Nuclear Star Clusters (NSCs) or “resolved 
stellar nuclei”. Such systems are observed in galaxies 
along the whole Hubble sequence (Richstone et al. 1998; 

Boker et al. 2002; Cote et al. 2006; Turner et al. 2012; den 
Brok et al. 2014) and they are the densest stellar systems 
known in the Universe, due to their high masses (up to 
10®Mq) and small scale radii (only few pc). Moreover, 
it is not rare to find in galactic nuclei both a NSC and a 
super massive black hole (SMBH) (Seth et al. 2008; Gra¬ 
ham & Spitler 2009). In particular, the center of galaxies 
with masses of ^ 10®Mq are typically dominated by the 
presence of a NSC, while a MBH could be present but it 
is likely unseen due to resolution limit of present instru¬ 
ments. On the other hand, galaxies with masses above 
1O**M0 are dominated by the presence of SMBHs. Be¬ 
tween such ranges, instead, SMBHs and NSCs have simi¬ 
lar masses (Graham & Spitler 2009; Neumayer & Walcher 
2012; Kormendy & Ho 2013). 

The fact that NSCs seem to dominate the center of 
smaller galaxies while SMBHs dominate heavier galac- 

m.arcasedda@gmail.com 


tic centers suggests a link between these objects which, 
however, remains poorly understood. 

In the last years, much research has been devoted to ex¬ 
amining scaling relations connecting the SMBHs or NSCs 
and their hosts. For instance, Ferrarese et al. (2006) and 
Rossa et al. (2006) have shown that the relation between 
the NSC mass and the host velocity dispersion is similar 
to that provided for SMBHs. On the other hand, more 
recent studies have shown that this relation is shallower 
for NSCs (Leigh et al. 2012; Erwin & Gadotti 2012; Scott 
& Graham 2013). The correlation between the central 
object (SMBH or NSC) and the host could provide im¬ 
portant information about the central object’s formation 
history and evolution. 

At present there are several formation models for NSCs. 
In the scenario commonly referred to as “in-situ” forma¬ 
tion (King 2003, 2005; Milosavljevic 2004; Bekki et al. 
2006; Aharon & Perets 2015), the NSC forms if the time- 
scale on which the BH can grow by accretion is larger 
than the crossing time of the galaxy (Nayakshin et al. 
2009; Hopkins & Quataert 2010). Such a scenario pro¬ 
vides a relation between the NSC mass and the host ve¬ 
locity dispersion similar to that observed for SMBHs. 
Other works proposed that a NSC forms first in the galac¬ 
tic nucleus, and then its high density facilitate runaway 
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collisions of massive stars, leading to the formation of a 
SMBH seed, which may subsequently grow through stel¬ 
lar and gas accretion (Portegies Zwart et al. 2004). 

On the other hand, NSCs can form via the well- 
studied “dry-merger” scenario, in which star clusters or- 
bitally decay toward the galactic center under the ac¬ 
tion of the dynamical friction process (Tremaine et al. 
1975; Tremaine 1976; Capuzzo-Dolcetta 1993; Capuzzo- 
Dolcetta & Mastrobuono-Battisti 2009). The subsequent 
merging of the decayed star clusters leads to the forma¬ 
tion of a dense stellar nucleus with characteristics com¬ 
parable to those of real nuclei. In this scenario, the pres¬ 
ence of a central BH in the host galaxy could lead to 
the disruption of the decaying clusters, preventing NSC 
formation. However, Antonini (2013) and Arca-Sedda 
& Capuzzo-Dolcetta (2014b) have shown that the tidal 
forces induced by the central SMBH are only significant 
at BH masses above IO^Mq. 

The dry-merger scenario has been investigated through 
numerical simulations (Capuzzo-Dolcetta & Miocchi 
2008b, a; Antonini et al. 2012; Perets & Mastrobuono- 
Battisti 2014), and also via statistical and theoretical 
models (Leigh et al. 2012; Antonini 2013; Arca-Sedda 
& Capuzzo-Dolcetta 2014b). These works have shown 
that such a scenario could explain the observed proper¬ 
ties of the nuclei, and also provide an explanation for 
the observed NSC-host galaxy scaling relations recently 
updated by several works (Leigh et al. 2012; Erwin & 
Gadotti 2012; Scott & Graham 2013). 

In this context, galaxies with star clusters near their pho¬ 
tometric center provide a perfect laboratory to test the 
dry-merger scenario. Two examples of such galaxies are 
the Fornax dSph galaxy, the heaviest satellite galaxy of 
the Milky Way (Buonanno et al. 1999, 1998; Mackey 
& Gilmore 2003; Dinescu et al. 2004; Goleman et al. 
2005; Goerdt et al. 2006; Read et al. 2006; Cole et al. 
2012), and the Henize 2-10 dwarf starburst galaxy (Kob- 
ulnicky et al. 1995; Johnson et al. 2000; Johnson & Kob- 
ulnicky 2003; Ott et al. 2005a, b; Santangelo et al. 2009; 
Reines & Deller 2012). A recent paper provided detailed 
constraints on the galaxy and star cluster properties of 
Henize 2-10 (Nguyen et al. 2014). This galaxy likely hosts 
an accreting SMBH (Merloni et al. 2003; Kobulnicky & 
Martin 2010) of mass Log{M/M q) = 6.3 ±1.1 (Reines 
et al. 2011) and eleven young super star clusters (SSCs) 
with masses above 10® Mq, placed at projected distances 
< 140 pc from the galactic center (Nguyen et al. 2014). 
On the other hand, no NSC has been detected. The de¬ 
tailed informations available on the projected mass pro¬ 
file of the galaxy and the cluster masses and sizes, make 
Henize 2-10 an interesting case study for testing the dry- 
merger scenario. 

In this paper we make use of direct-summation 7V-body 
simulations to investigate whether the future evolution 
of the young star cluster system observed in Henize 2-10 
could lead to the formation of a central stellar nucleus, 
distributing the SSCs with initial conditions consistent 
with observational constraints on their orbit and inter¬ 
nal properties. 

The paper is organized as follows: in Section 2 we in¬ 
troduce the models used to represent galaxy and SSCs 
and give an estimate of the mass growth of the expected 
NSC by using proper semi-analytical arguments; in Sec- 


Table 1 

Adopted parameters of the super-star clusters. 


ID 

Mssc 

(10®Mq) 

rc 

(pc) 

reff 

(pc) 

Rssc 

(pc) 

N 

G1 

2.30 ±0.60 

1.27 ±0.03 

3.1 

64.8 ± 1.4 

28976 

C2 

0.92 ± 0.61 

0.77 ±0.01 

1.9 

36.9 ± 1.4 

11590 

C3 

1.14 ±0.78 

1.00 ±0.02 

2.4 

20.4 ±1.4 

14362 

C4 

0.91 ± 0.60 

0.71 ±0.05 

4.2 

11.2 ± 1.4 

11464 

C5 

0.40 ±0.21 

0.53 ±0.07 

1.3 

86.4 ±1.4 

5039 

C6 

0.40 ±0.21 

0.53 ±0.04 

1.3 

77.6 ± 1.4 

5039 

C7 

0.46 ±0.25 

0.89 ±0.21 

2.2 

114.5 ± 1.4 

5795 

C8 

0.45 ±0.24 

0.46 ±0.02 

1.1 

130.7 ±1.4 

5669 

C9 

0.20 ±0.11 

0.37 ±0.01 

0.9 

74.1 ± 1.4 

2519 

CIO 

0.45 ±0.24 

0.46 ±0.02 

1.1 

70.2 ± 1.4 

5669 

Cll 

0.20 ±0.11 

0.37 ±0.01 

0.9 

55.8 ± 1.4 

2519 


Column 1: cluster name. Column 2: mass. Column 3: core radius. 
Column 4: effective radius. Column 5: projected radial coordinate. 
Column 6: number of particles used. 

tion 3 the results of the simulations are presented and 
discussed; in Section 4 we discuss the possible long-term 
evolution of the NSC in Henize 2-10; finally. Section 5 is 
devoted to a final discussion and conclusions. 

2. MODELING THE GALAXY AND THE GLOBULAR 
GLUSTER SYSTEM 

2.1. The galaxy model 

Data provided by Nguyen et al. (2014) give virial 
estimates of the mass enclosed within two different 
values of the distance to the center of Henize 2-10: 

M(r < 2”) = (2.7 ± 1.1) X 10®Mq, (1) 

M(r < 6”) = (6.4 ± 1.5) X 10®Mq; (2) 

as well as the effective radius of the inner component of 
the galaxy, Re = 6”. At the estimated distance of 9.0 
Mpc to us (Vacca & Conti 1992), the above values refer 
to the mass enclosed within 89 and 259 pc, respectively. 
We modeled the galaxy using the density profile model: 

(3-7)Mg (r 1 

^ 47rr® Vrs/ \rs ) cosh(r/rtr) ’ 

( 3 ) 

where Mh is the total mass of the galaxy, Vg its scale 
radius, 7 is the inner density profile slope, and rtr is the 
truncation radius of the model. At radii r ^ rtj. the 
model corresponds to a Dehnen density profile (Dehnen 
1993) with mass distribution: 

M,(r) = MH(^) ", (4) 

The scale length, r^, is related to the observed effective 
radius, i?e, through the relation 

^(2^/(3-^) -l)Re. (5) 


To limit the number of particles in our direct 7V-body 
simulations of Henize 2-10 we decided to restrict our self- 
consistent, particle, representation of the galaxy to the 
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Figure 1. Cumulative projected mass distribution (solid line) as 
obtained by Equation 3 with the values given in Subsect. 2.1, com¬ 
pared with the measured values of the mass enclosed within 2” 
(86 pc) and 6” (259 pc). The vertical dashed line represents the 
effective radius of the galaxy. 

inner radial region. We show that this simplification pro¬ 
vides a good description of the environment in which the 
clusters move and discuss the resolution of the simula¬ 
tions in Section 3.1. 

Due to this choice, in what follows we set Mh = 1.6 x 10® 
Mq, 7 = 1/2 and rtr = 150 pc. This choice accu¬ 
rately represents the deprojected brightness profile of 
the galaxy at r < 300 pc and leads to a value of the 
galaxy mass within 6 ” in good agreement with Equa¬ 
tion 2. Moreover, the projected velocity dispersion in 
the range 1” — 6 ” has a value of 41 km s“^, very close to 
the observed value (Marquart et al. 2007; Nguyen et al. 
2014), ensuring a reliable representation of the dynam¬ 
ics of the clusters within the galaxy. It is worth noting 
that the high resolution of the data provided in Nguyen 
et al. (2014), ^ 0.1” , implies that an error on the po¬ 
sitional center of the galaxy does not affect significantly 
our choice of the galaxy model. 

In Figure 1 we compare the mass profile of our model 
with the mass profile obtained from the photometric data 
of Nguyen et al. (2014), assuming a distance of d = 9.0 
Mpc (Vacca & Conti 1992). 

Since there is observational evidence that the galaxy 
hosts a MBH at its center with mass Log(MBH/M 0 ) = 
6.3±1.1, we decided to include it in our models, generat¬ 
ing the A^-body representation of the whole galaxy via a 
numerically calculated distribution function. While the 
mass of the BH candidate is quite uncertain, and there¬ 
fore it would be worth testing models with different BH 
masses, the numerical simulations presented here are ex¬ 
tremely time consuming. Therefore, we examine only one 
value for the BH mass (Mbh = 2.6 x 10® M 0 ), while vary¬ 
ing the initial conditions for the clusters to explore the 
future evolution of the galactic nucleus. Furthermore, we 
ran also one simulation in which the galaxy does not host 
any MBH, in order to highlight how a central, massive 
object can affect the formation of a bright nucleus. 

2 . 2 . The globular cluster system model 

We modeled each of the Henize 2-10 SSCs using a 
King model (King 1966), for which Nguyen et al. (2014) 
provides mass, Mgsc, core radius, rc, and projected 
galactocentric distance, Rssc- We investigate three 
different possibilities for the initial orbital distribution 


of the 11 SSCs. In the first two models we vary initial 
conditions for the SSCs, while, in the third case, clusters 
are distributed as in the first simulation, but the galaxy 
does not contain a central BH. The first two initial 
conditions provide bounds on the dynamical friction 
inspiral times of the clusters. In what follows we will 
refer to these three choices of initial conditions as models 
SI, S2 and S3. In model SI the clusters were assumed 
to follow the distribution of the background galaxy, but 
with the constraint that the clusters projected position 
(onto a random plane) was consistent with the observed 
values of Rssc- In order to do so we first sample the 
galaxy using the density model of Equation 3, and 
then selected points from that sample with projected 
positions (on a random plane) similar to those observed 
for Henize 2-10 SSCs. The clusters initial velocities 
were drawn from an isotropic distribution function 
corresponding to the density model of Equation 3. This 
model represents the one in which the action of the 
dynamical friction is minimized, since the clusters have 
initial positions that exceed their projected (observed) 
positions. 

Model S2 corresponds to the assumption that the 
clusters are all located on the same plane, so that the 
projected and the spatial distances coincide, and that 
they are on circular orbits with velocities as required 
by the mass distribution of the galaxy. This scenario 
is motivated by the observed velocities of the clusters 
that rotate in the same sense as the molecular gas in 
Henize 2-10, and thus may lie in a disk (Santangelo 
et al. 2009; Nguyen et al. 2014). In this case the 
projected and spatial positions of the clusters coincide, 
and therefore the effects due to the dynamical friction 
are the strongest possible. 

Finally, for the model S3 we used the same initial 
conditions as in SI model, but in this case the galaxy 
does not contain a central BH. This simulation allows us 
to study the effect of an SMBH on the NSC formation 
process. 

Once the spatial positions of the clusters were obtained, 
we evaluated their tidal radius, r^, a parameter needed 
to define the cluster King model. In the hypothesis 
of circular orbit, the tidal radius of the cluster is given by: 

^3 ^_ GMssc _ 

* + (d2t//dr2),33, ’ ^ ^ 

where w is the angular velocity of the cluster and 
{d^U/dr^)rssc is the second derivative of the potential 
evaluated along the SSC circular orbit (rssc = const.). 
The knowledge of rt and Tc makes possible the evaluation 
of the concentration of the cluster, c = rtjrc^ which cor¬ 
relates with the dimensionless potential well depth, Wq, 
needed to dehne the King model. 

Another parameter we can compare with observations is 
the effective radius of the cluster, reff, which can be eval¬ 
uated from Tc through the relation (Miocchi et al. 2013): 

reff = (2.43 ± 0.9)re. (7) 

The effective radii for the modeled cluster evaluated this 
way and reported in Table 1 are in agreement with those 
in Nguyen et al. (2014). 

Tables 2 and 3 provide the main parameters for the star 
cluster system (SCS) models, SI, S2, and S3 considered 
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Figure 2. Mass accreted around the galactic center as function 
of time assuming the dynamical friction formula presented in Eq. 
8. Different colors correspond to different set of initial conditions 
for the clusters. The NSC is expected to form within ~ 500 Myr. 

Table 2 

Main parameters of SSCs in configurations SI and S3. 


Table 3 

Main parameters of SSCs in configuration S2. 


ID 

Wo 

c 

n 

(pc) 

(7 

(km s“^) 

rsSC 

(pc) 

Ta 

(pc) 

e 

T-df 

(Myr) 

Cl 

5.1 

11.5 

14.6 

26.1 

64.3 

64.3 

0 

20.0 

C2 

5.1 

11.6 

8.9 

21.1 

37.0 

37.0 

0 

13.9 

C3 

4.2 

7.6 

7.6 

25.3 

20.5 

20.5 

0 

4.3 

C4 

4.0 

6.9 

4.9 

28.2 

11.3 

11.3 

0 

1.7 

C5 

5.9 

17.1 

9.1 

13.8 

85.6 

85.7 

0 

106.7 

C6 

5.8 

16.5 

8.7 

14.0 

78.0 

78.0 

0 

90.6 

C7 

5.2 

12.2 

10.8 

13.5 

114.2 

114.2 

0 

161.3 

C8 

6.4 

25.0 

11.5 

13.0 

131.3 

131.3 

0 

209.3 

C9 

6.0 

18.4 

6.8 

11.2 

74.5 

74.5 

0 

132.9 

CIO 

6.1 

19.0 

00 

bo 

14.9 

70.9 

70.9 

0 

70.7 

Cll 

5.8 

16.6 

6.1 

11.8 

56.1 

56.1 

0 

80.7 


Column 1: cluster name. Column 2: adimensional central poten¬ 
tial. Column 3: concentration parameter. Column 4: tidal radius. 
Column 5: velocity dispersion. Column 6: galactocentric distance. 
Column 7: initial apocentric distance of the SSC orbit. Column 
8: eccentricity of the orbit. Column 9: dynamical friction time 
evaluated using Eq. 8. 


ID 

Wo 

c 

rt 

(pc) 

CT 

(km s“^) 

rsSC 

(pc) 

ra 

(pc) 

e 

(Myr) 

Cl 

5.3 

12.9 

16.5 

24.5 

87.9 

145 

0.4 

43.4 

C2 

5.5 

14.5 

11.2 

18.8 

71.6 

75 

0.3 

28.4 

C3 

5.1 

11.7 

11.7 

20.5 

67.6 

200 

0.8 

85.3 

C4 

4.1 

7.4 

5.3 

27.3 

12.4 

30 

0.4 

7.1 

C5 

6.1 

19.1 

10.1 

13.1 

109.1 

120 

0.5 

116.2 

C6 

6.2 

22.3 

11.8 

12.1 

148.1 

190 

0.5 

270.0 

C7 

5.2 

12.3 

10.9 

13.5 

116.4 

140 

0.5 

138.8 

C8 

6.5 

25.2 

11.6 

12.9 

132.5 

160 

0.6 

163.4 

C9 

6.2 

21.0 

7.8 

10.5 

101.6 

180 

0.5 

377.4 

CIO 

6.1 

19.2 

8.8 

14.8 

72.2 

180 

0.6 

203.0 

Cll 

5.9 

17.7 

6.5 

11.5 

66.6 

140 

0.4 

265.4 


Column 1: cluster name. Column 2: adimensional central poten¬ 
tial. Column 3: concentration parameter. Column 4: tidal radius. 
Column 5: velocity dispersion. Column 6: galactocentric distance. 
Column 7: initial apocentric distance of the SSC orbit. Column 
8: eccentricity of the orbit. Column 9: dynamical friction time 
evaluated using Eq. 8. 

in this paper. Note that the cluster parameters for sim- 
ulations SI and S3, are the same since the positions of 
the clusters are the same in these simulations. 

2.3. The NSC formation process in Henize 2-10: 
analytical predictions 

Dynamical friction leads star clusters to decay toward 
the center of their parent galaxies (see for example 
Tremaine et al. (1975), Capuzzo-Dolcetta (1993)). 
In general, the easiest way to describe the action of 
dynamical friction in a galactic environment is by means 
of Chandrasekhar’s formula (Chandrasekhar 1943) in its 
local approximation. However, many works (Ostriker 
et al. 1989; Pesce et al. 1992; Capuzzo-Dolcetta & 
Vicari 2005; Arca-Sedda & Capuzzo-Dolcetta 2014a) 
have shown that this local approximation fails when the 
satellite moves in the innermost region of a galaxy. 
Arca-Sedda & Capuzzo-Dolcetta (2014a) developed a 
useful formula for the dynamical friction decay time 
valid for both cored and cuspy density profiles, and for 
galaxies with a central black hole: 


rdf(Myr) = 0.3 



/ Mh \ 
\Mssg) 



( 8 ) 


where 7 is the previously defined exponent in the galac¬ 
tic density profile, 0 < e < 1 is the orbital eccentricity 
of the point-like object whose mass is Mssc, fa. is the 
orbit apocenter. Here, we used for the eccentricity the 
definition: 



being ra and rp the apocentral and pericentral dis¬ 
tance, respectively. Parameters rg (in kpc) and Mh (in 
IO^^Mq) are the scale radius and the total mass of the 
galaxy 7 model. For the scopes of this paper, we per¬ 
formed several new direct N-body runs to test the va¬ 
lidity of the formula for massive objects on circular and 
quasi radial orbits and for 7 values smaller than 2 , find¬ 
ing that a function 5 ( 6 , 7 ) that gives a very good ht to 
results has this expression: 


5 ( 6 , 7 ) = (2 - 7 ) 


Ol 


(2 - 7 )“ 


-I- 03 I (1 - 6) -I- 6 


( 10 ) 


with Ol = 2.63 ± 0.17, 02 = 2.26 ± 0.08 and 03 = 
0.9 ± 0.1. The expression for 5 ( 6 , 7 ) here is slightly dif¬ 
ferent from the one found in Arca-Sedda & Capuzzo- 
Dolcetta (2014a) due to that that paper the authors 
mainly focused their attention on the study of the dy¬ 
namical friction process in models with 7 < 1 for both 
circular and radial orbits, and 7 > 1 limiting to nearly 
radial orbits. Here, we have extended the interval of 
validity of the fitting formula to circular and radial sim¬ 
ulations in the range 0 < 7 < 2. Note, however, that 
Equation 8 is valid only for 7 < 2. 

Equation 8 can be used to estimate the mass accumula¬ 
tion around the center of Henize 2-10 as a consequence 
of the orbital decay of the 11 clusters. 

To have a meaningful idea of what we should expect from 
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Table 4 

Estimated initial galactocentric distances of the SSCs for different 
assumptions of their ages. 


ID 

roAssC 


Cl 

1.1 

1.3 

2.8 

C2 

1.2 

1.4 

3.3 

C3 

1.6 

2.0 

6.1 

C4 

2.2 

3.0 

10.1 

C5 

1.0 

1.1 

1.5 

C6 

1.0 

1.1 

1.5 

C7 

1.0 

1.0 

1.3 

C8 

1.0 

1.0 

1.2 

C9 

1.0 

1.0 

1.4 

CIO 

1.0 

1.1 

1.7 

Cll 

1.0 

1.1 

1.6 


Column 1: cluster name. Column 2-4: ratio between the initial 
galactocentric distance the cluster had at its birth, ro, and its 
current position, rssCi assuming as age estimates: 5, 10 and 100 
Myr, respectively. 

iV-body simulations, for each of the 11 SSCs in Henize 
2-10 we randomly selected initial orbital eccentricity and 
position values. Then, for any given set of initial parame¬ 
ters we evaluated the mass accreted to the galaxy center. 
Figure 2 shows the mass accumulated around the center 
of the galaxy for 7 sets of random initial conditions: a 
mass of ^ 7x10® M© is accreted into the center of the 
galaxy with timescales between 0.1 and 0.5 Gyr. This 
mass is consistent with the mass of observed NSCs in 
other galaxies (e.g. Seth et al. 2008). 

Tables 2 and 3 give the decay time of each cluster in 
the configuration SI and S2, respectively, evaluated us¬ 
ing Equation 8. 

The semi-analytical formula of Equation 8 neglects the 
mass loss of clusters along their orbits, and thus repre¬ 
sents a lower limit to the time needed to the decay and 
an upper limit to the quantity of cluster mass deposited 
into a NSC. A more realistic treatment requires N-body 
simulations, which we present in the following sections. 

The work by Chandar et al. (2003) suggests that all 
the 11 SSCs studied are very young, r ~ 5 Myr, therefore 
their initial galactocentric distances, rg, are likely very 
close to the presently observed, rssc, values. We can 
quantify this by mean of an easy inversion of Equation 
8 , which tells us how much the SSCs orbits have shrunk 
over their (short) life. The data Table 4 are obtained 
in the simplifying assumption of SSCs on initially circu¬ 
lar orbits with following circular orbital shrink. Data in 
Table 4 indicate, even in the very unlikely case that the 
clusters age reaches 100 Myr, they should have formed 
within 1 kpc from the center of Henize 2-10. 

3. RESULTS FROM W-BODY SIMULATIONS 
3.1. The “number of particles” problem 

In this work we simulate the future dynamical evo¬ 
lution of the clusters observed in Henize 2-10, to see 
whether or not a nuclear stellar cluster will form, and 
what the timescale of this formation process is. To do 
this, we used the direct V-body code HiCPUs (Capuzzo- 
Dolcetta et al. 2013), a 6‘^-order Hermite integrator that 
runs on hybrid platforms where CPUs are coupled to 
Graphic Processing Units (GPUs) acting as computing 
accelerators, fully exploiting parallelization. 


Despite the great accuracy and speed of the code, the 
number of particles that we can use in this kind of simula¬ 
tions is limited to ^ 10®. Since Henize 2-10 has a mass of 
10 ^® Mq, it is clearly impossible to create a reliable, star- 
by-star simulation of the whole galaxy. To address this 
problem, we limit the extension of our iV-body sampling 
of the galaxy to a relatively small region, as explained 
in Section 2.1. In the simulations we set a gravitational 
softening e = 0.02 pc, which is a value much smaller than 
the core radius of the clusters, and thus allows for a good 
description of the internal dynamics of the clusters. 

Moreover, we allowed for a difference between the mass 
of particles representing the galaxy, mg, and that of the 
particles in the clusters, mssCi such that mg/mssc = 8. 
This allows us to have a number of particles in the clus¬ 
ters, with even the smallest cluster having N ~ 2000, 
sufficient to ensure reliable results on their dynamical 
evolution. Specifically, the two-body relaxation time of 
the smallest cluster is long enough to make the evapora¬ 
tion time Tevp — IdOTrei — 2.5 Gyr (Binney & Tremaine 
2008), much longer than the run time of our Wbody 
simulations. This guarantees that the global properties 
of the small clusters are preserved along the integration 
time despite being undersampled. 

3.2. Configuration SI 

As explained above, the clusters in the SI configuration 
have initial conditions obtained directly from the distri¬ 
bution function of the galaxy model given in Equation 
3. 

Table 2 summarizes the parameters that describe the 
SCS in this case. 

In the following, each cluster will be identified with the 
letter C and a number between 1 and 11, as in Table 1. 

Figure 3 shows the distance of each cluster from the 
central galactic BH as a function of time, as obtained by 
our iV-body simulation. It is worth noting that an error 
on the position of the galaxy center would not affect our 
estimates, since we are interested in the mutual distances 
between the clusters and the BH. As expected, while the 
most massive clusters (C1-C4) reach the BH within ^ 50 
Myr, the dynamical friction is less efficient on lighter 
clusters, which tend to remain on their initial orbits for 
longer times. 

Figure 4 shows that the time evolution of the cluster or¬ 
bits obtained through the A^-body simulation agrees with 
the semi-analytical method described in Arca-Sedda & 
Capuzzo-Dolcetta (2014a). It is evident that the semi- 
analytical orbit accurately describe the cluster orbital 
evolution for the most massive clusters, confirming the 
reliability of Equation 8. On the other hand, the motion 
of lighter clusters is strongly affected by the tidal forces 
induced by the BH and the other clusters. Therefore, in 
this case the semi-analytical approach is less accurate, 
but even for these clusters the results are in broad agree¬ 
ment with the 7V-body results. 

As the clusters move within the galaxy, their initial 
spherical shape changes significantly due to the action of 
the tidal forces. In this case, it is very difficult to define 
general quantities of the clusters such as their center or 
mass. To give an estimate of the cluster mass at any time, 
we obtained rt solving Equation 6, while the position 
of the cluster center is estimated through a grid-based 
algorithm. Once the tidal radius and cluster center have 
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Figure 3. Left panels: radial distance of each cluster from the BH as a function of time. Right panels: mass of each cluster as function of 
time. It is evident that the main contributors to the nuclear mass are clusters C1-C4, which transport a significant amount of their initial 
mass to the center. Panels refer to simulation SI 


been evaluated, we considered the cluster mass in a given 
time-step to be the mass enclosed within 

In the right panels of Figure 3 the evolution in time of 
the mass of the clusters is shown. In this case, it is evi¬ 
dent that massive clusters (C1-C4), lose their mass when 
they reach the galaxy center, contributing signihcantly 
to the formation of the nucleus. 

On the other hand, small clusters (C5-C10) lose mass 
continuously throughout the simulation to the field. The 
evolution of cluster C8 and Cll are significantly affected 
by tidal forces, which strip away more than 70% of the 
initial clusters mass well before they reach the center of 
the galaxy. 

Combining the informations on the time evolution of 


the SSCs mass and trajectory, we can extrapolate the 
amount of mass that each cluster will carry to the inner¬ 
most region of the galaxy. In particular, Table 5 shows 
the percentage of mass dragged within 20 pc from the 
galactic center for all the simulations performed. In the 
case of configuration SI, it is evident that tidal forces 
should prevent the complete decay of small clusters, dis¬ 
rupting completely clusters C5-C9 and Cll, and leaving 
cluster CIO with only ^ 30% of its initial mass. 

Figure 5 gives information on the mass transported to 
the galactic center, showing the cumulative radial mass 
distribution of the SCS at different times. It is evident an 
initial phase lasting ^ 40 Myr in which the mass accumu¬ 
lated within 20 pc increases until it reaches 55% of the 
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Figure 4. From top to bottom: radial distance from the BH of 
the cluster Cl, C3 and C5 as a function of time obtained with 
the self-consistent simulation SI (thick black line), and the semi- 
analytical method described in Arca-Sedda &; Capuzzo-Dolcetta 
(2014a) (thin red line). The comparison shows a good agreement 
for heavy clusters, while the analytical description is in poorer 
agreement for smaller clusters. 

total SCS mass. This represents in this conhguration a 
“saturation” value for the accumulated mass. Combining 
this figure with Figure 3 reveals that the major contri¬ 
bution to the central mass is from the C1-C4 clusters, 
while the lighter clusters decay more slowly. Moreover, 
tidal forces act more efficiently on small clusters, leading 
to mass loss; this makes the decay of these clusters even 
slower. 

Figure 6 shows the variation with time of the distance 
of the BH to the galactic center, revealing that in this 
SI case the BH motion remains limited to a radial re¬ 
gion of <6 pc from its initial position. The maximum 


Table 5 

Percentage of mass left to the galactic center by each cluster in 
the three configurations studied. 


ID 

SI 

S2 

S3 

Cl 

80 

80 

93 

C2 

80 

95 

89 

C3 

70 

95 

73 

C4 

95 

95 

98 

C5 

0 

30 

0 

C6 

0 

40 

0 

C7 

0 

40 

0 

C8 

0 

10 

0 

C9 

0 

0 

0 

CIO 

30 

70 

30 

Cll 

0 

60 

0 


Column 1: cluster name. Column 2-4: percentage of the cluster 
mass left within 20 pc from the BH in configuration SI, S2 and S3, 
respectively. 



Figure 5. The mass distribution, as function of the projected 
radial distance from the BH, of the SCS at different times, normal¬ 
ized to the total SCS mass for configuration SI. The mass enclosed 
within 20 pc reaches a nearly constant value after 40 Myr. 



Figure 6 . Evolution of the distance of the BH from the galaxy 
center as function of time in cases SI and S2. 

displacement, at about 40 Myr, corresponds to the close 
interaction with Cl (the heaviest cluster), as shown in a 
comparison of Figure 6 and Figure 3. After the comple¬ 
tion of the main mergers, the BH oscillates within 1 pc 
from the center with a residual speed of ^ 2 km s“^ (see 
Figure 7). 
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Figure 7. Time evolution of the BH speed in cases SI and S2. 



Figure 8 . Mass deposited in a sphere centered in the BH of radius: 
4 pc (straight line), 10 pc (dashed line), and 20 pc (dotted line), 
for the simulation SI. 

To understand whether a NSC forms, we evaluated 
the cluster mass deposited in the innermost region of 
the galaxy. This is shown in Figure 8, which shows the 
amount of cluster mass within 4, 10 and 20 pc distance 
from the central BH as a function of time. 

The evident steps in the increasing form of the function 
flag the time at which a new cluster reaches the Henize 
center. After 90 Myr, the deposited mass within 20 pc 
reaches a nearly constant value of 4 x 10® Mq, i.e. 51% 
of the total SCS mass. 

The presence of a NSC in the center of a galaxy is, ob- 
servationally, usually revealed by the study of the surface 
luminosity profile of the inner part of the galaxy. In fact, 
a NSC tipically emerges in such profiles as an overdensity 
in the innermost region. 

The initial and final surface density profiles are shown 
in Figure 9. It shows clearly the formation of a well- 
visible structure which extends up to 0.2” (10 pc). Since 
the detection of a NSC is made observationally by eval¬ 
uating the mass contained within the overdensity, our 
newly born NSC has an observational extension of 10 
pc, a mass of Mnsc = 4.6 x 10® Mq and an effective ra¬ 
dius tnsc — 4.17 pc. It should be highlighted that, here, 
Mnsc is obtained summing the mass which the decayed 
clusters carried to the galactic center and the mass of the 
galactic background enclosed within 10 pc from the BH. 
On the other hand, the mass deposited through the clus¬ 
ter decay in the range 10 — 20 pc (< 2 x 10® Mq) is 



Figure 9. Surface density profile of the system galaxy+SCS at 
time t = 0 Myr (straight line) and t = 83.1 Myr (dotted line) for 
the configuration SI. The x-axis is the projected distance from the 
galaxy center. 

evidently not dense enough to give rise to a clearly de¬ 
tectable structure. 

3.3. Configuration S2 

In the conhguration S2 the SSC orbits lie on the plane 
perpendicular to the line of sight. Such a configuration 
would result in the strongest dynamical friction. This 
scenario may also be closest to the clusters true distribu¬ 
tion. Indeed, the young clusters, all located to the East 
of the black hole and dynamical center, are all blueshifted 
(Nguyen et al. 2014). The gas in this region is rotating 
in the same direction (Cresci et al. 2010), and overall the 
gas in the galaxy has a low inclination, i ^38° (Kobul- 
nicky et al. 1995). Thus it appears reasonable that the 
young clusters are in a rotating disk close to the plane of 
the sky. 

Parameters of the clusters in this configuration are 
summarized in Table 3. 

In this case, clusters reach the galactic center in shorter 
times than in case SI since they move on projected orbits, 
which are of course smaller than the corresponding 3D 
positions. 

The exceptions are clusters C7 and C8, that move as 
a binary. 

Moreover, clusters that reach the center transport a 
significant fraction of their mass toward the nucleus, as 
shown in the right panels of Figure 10 and in Table 5. 

The distance of the BH from the galactic center as a 
function of time is shown in Figure 6. After some initial 
oscillations around its initial position with amplitude ^ 
2 pc, the BH returns to the galactic center within 20 
Myr. The damping of the oscillations is due mainly to 
dynamical friction induced by the background stars. 

The cumulative radial mass distribution of the SCS 
at different times in this S2 model is reported in Figure 
11. It is interesting to note that after ^ 30 Myr the 
mass distribution enclosed within 20 pc from the BH 
does not change significantly since the dynamical friction 
timescale in this configuration does not exceed 20 Myr 
for all the clusters. The only exception is cluster Cll, 
whose decay time is 40 Myr. 

The mass deposited within 4, 10 and 20 pc from the 
BH is shown in Figure 12. The Figure makes evident 
that, in this case, more than 50% of the SCS total mass 
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Figure 10. As in Figure 3, but for simulation S2. In this case all 

reaches the galactic center giving rise to a stellar nucleus, 
clearly visible out to 10 pc (0.2”) in the surface density 
profile, as shown in Figure 13. 

Such a structure has a total mass Mnsc(?' < 10 pc) ~ 
6 X 10® Mq and an effective radius msc — 2.63 pc. 

3.4. Configuration S3 

In the simulations presented above we studied the fu¬ 
ture evolution of the Henize 2-10 central region, assuming 
the presence of a central BH with mass Mbh = 2.6 x 10® 
Mq. However, the mass measurements from Heines et al. 
(2011) are based on the fundamental plane and thus quite 
uncertain, placing a lower limit on the mass of the BH 
candidate at Mbh = 1.6 x 10® Mq. More recent results 
by Nguyen et al. (2014) give an upper limit of the BH 





0 10 20 30 40 50 60 70 80 

t (Myr) 

the clusters contribute to the accretion of a nucleus. 

mass of Mbh = lO"^ M©. Therefore, it is relevant under¬ 
standing the role of variation of the black hole mass on 
the SSC evolution. 

As pointed out by several authors (Capuzzo-Dolcetta 
1993; Antonini 2013; Arca-Sedda & Capuzzo-Dolcetta 
2014a; Spera et al. 2015), the main effects of a central 
BH on the motion of infalling star clusters are that of 
a slight reduction of the dynamical friction decay time, 
enhancing at the same time the tidal disruption of the 
infalling clusters. 

As shown in Arca-Sedda & Capuzzo-Dolcetta (2014a), 
a central BH may reduce the decay time if: i) the 
orbit of the cluster is not far from radial, and ii) if 
Mssc/Mqh ~ 1. Otherwise, the presence of a BH does 
not alter significantly the clusters decay times. 
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Figure 11. Cumulative projected mass distribution of the SCS 
at different times in the configuration S2. The x-axis reports the 
projected distance from the galactic center. Here the deposited 
mass reach a saturation value on time-scales smaller than 20 Myr. 



Figure 12. The same as in Figure 8, but for the configuration S2. 



Figure 13. Initial (straight line) and final (dotted line) surface 
density profile of the system galaxyH-SCS. The x-axis reports the 
projected distance from the galactic center. In this case the nucleus 
emerges from the galactic background. 

On the other hand, the tidal forces induced by the BH 
may quench the NSC formation by cluster mergers, dis¬ 
rupting the incoming clusters before they reach the galac¬ 
tic center. Recent work has demonstrated that whenever 
the infalling cluster has a mass Mssc < O.OIMbh, the 
tidal forces induced by the BH efficiently disrupt the clus¬ 
ter as it moves within the BH influence radius, thus sup¬ 


pressing the formation of a NSC (Antonini 2013; Arca- 
Sedda & Capuzzo-Dolcetta 2014b; Spera et al. 2015). 

To investigate the quantitative role of the BH in the 
NSC formation process, we decided to study the extreme 
case of absence of a central BH in a numerical simulation 
where all the other initial conditions are as in model SI. 

The main effect of the actual presence of a MBH can be 
argued from Figure 14, which shows the time evolution of 
the radial distance of SSCs from the galactic center and 
of the SSCs masses. A comparison of this figure with 
Figure 3 shows how the MBH, on one side, affects just 
slightly the motion of the SSCs, but, on another side, it 
plays a significant role in determining the amount of mass 
dragged to the center of the host galaxy. Indeed, while 
the time-scale for the orbital decay is almost the same 
in the two (SI and S3) configurations, we found that in 
model S3 the SSCs reach the galactic center keeping a 
mass larger than in SI. For instance, in S3 clusters Cl 
and C2 reach the Henize 2-10 center conserving, respec¬ 
tively, more than the 93% and 89% of their initial masses, 
while in simulation SI the percentage of the bound mass 
is 80% for cluster Cl and ^ 80% for cluster C2. This is 
likely due to the tidal forces induced by the MBH, which 
facilitate the mass erosion of the infalling clusters. 

In Table 5 we listed the fraction of mass of each cluster 
left around the galactic center in configuration S3. Com¬ 
paring such values with those of SI makes clear that the 
most massive clusters, which reach earlier the galaxy cen¬ 
ter, have final masses larger than those in configuration 
SI while, on the other hand, lighter clusters are almost 
completely disrupted before they get to the innermost 
region of the galaxy, with the only exception of cluster 
CIO, which keeps 30% of its initial mass by the end of 
the orbital decay in both the simulations. 

Such results can be interpreted as follows. The most mas¬ 
sive clusters, having shorter decay time, reach rapidly 
the galactic center, where their evolution is determined 
by the combined tidal action of the galactic background 
and of the MBH. On the other hand, lighter clusters 
spend most of their time relatively far from the galactic 
center, where the tidal forces induced by the MBH are 
less efficient. Due to this, the presence of an MBH does 
not influence relevantly the evolution of clusters C5-C9 
and Cll, which are progressively disrupted by the tidal 
forces exerted by the galactic background. This is fur¬ 
ther ly clarified by the evolution of cluster CIO, which 
reaches the galactic center with the same final mass in 
both simulation SI (where the MBH is present) and S3 
(where it is not). Actually, in the configuration SI the 
distance between CIO and the MBH never falls below 20 
pc, a region which encloses a stellar background mass 
(^ 1.5 X 10^ Mq) much larger than Mbh- Hence, the 
evolution of clusters with long decay times, which move 
sufficiently far from the MBH, is mostly determined by 
the tidal forces induced by the stellar background, while 
the effect due to the MBH is negligible. 

We try to resume and generalize our understanding 
of the effect of the presence and abscence of a MBH as 
studied in this paper. 

In our simulations, we showed that tidal heating leads 
to significant disruption of clusters C5-CII in all the 
cases studied. A BH candidate with a mass Mbh = 
10^ Mq will disrupt those clusters more efficiently, but 
clearly this does not change the amount of mass that 
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can be deposited toward the galactic center, since clus¬ 
ters C5-C11 are disrupted also by lighter BHs and only a 
very tiny fraction of their initial mass contributes to the 
growth of a bright nucleus. Moreover, in this case the 
heaviest clusters, C1-C4, have masses Mssc > O.OOMbh, 
and therefore tidal disruption process is not efficient 
enough to prevent the clusters decay, but may lead to a 
decrease in the accumulated matter of 10% (Spera et al. 
2015). 

As we showed above, simulation SI and S3 refer to the 
same set of initial conditions for the SCS, but SI model 
contains a central BH while S3 does not. The case in 
which Mbh = 1.6 X 10® Mq lies in this two possible 
models. Therefore, keeping models SI and S3 as terms 
of comparison, we expect that the amount of mass de¬ 
posited within 20 pc from the galactic center in such a 
case, should be in the range M^ep = 4 x 10® Mq (model 
SI) and Mdep = 4.5 x 10® Mq (model S3). 

Figure 15, which shows the cumulative projected mass 
distribution of the SCS at different times, provides fur¬ 
ther information about the amount of mass left to the 
galactic center. In particular, it is evident that, over a 
time-scale t > 40 Myr, the mass deposited within a pro¬ 
jected radius i? = 20 pc rises to a saturation value of 
^ 5 X 10® Mq. The subsequent evolution, instead, leads 
to a further concentration of the SCS. Indeed, at later 
times, most of the deposited mass moves toward the first 
10 pc of the galaxy. 

The mass deposited, Mnuc, within 4, 10 and 20 pc from 
the center of the galaxy is shown in Figure 16. A compar¬ 
ison with Figure 14 makes clear that the growth of Mnuc 
is due to the decay of the most massive clusters. Indeed, 
each step-like increase corresponds to the decay of one 
of the most massive clusters. In particular, the first in¬ 
crease at t ~ 7 Myr corresponds to the decay of cluster 
C4, the second to the decay of cluster C2, and, finally, 
the last two steps correspond to the decay of clusters Cl 
and C3, respectively. Furthermore, Figure 17 shows the 
initial and final surface density profile of the SCS. In this 
configuration, it is well evident the formation of a struc¬ 
ture that extends up to 6 pc (0.16”), enclosing a mass 
Mnsc — 5.1 X 10® Mq, and characterized by an effective 
radius t’nsc = 2 pc. 

4. HINTS ON THE POSSIBLE LONG-TERM EVOLUTION 
OF THE NUCLEUS 

Despite a detailed study of how the newly formed NSC 
and the surrounding nucleus will evolve and interact 
would require times too long to be simulated with high 
precision simulations, we can provide some hints about 
the subsequent evolution of the dense galactic center us¬ 
ing the results obtained above. 

In particular, in our iV-body simulations we found that 
the mass enclosed within 20 pc saturates to a value 
4 X 10® Mq for the simulation SI, ~ 5 x 10® Mq for the 
simulation S2 and ^ 5 x 10® Mq for the simulation S3. 

Moreover, it is worth noting that the clusters which 
still were on their orbit at the end of the simulations, are 
expected to give a negligible contribution to the growth 
of the nucleus, since they are likely disrupted before they 
can reach the inner region of the galaxy (see Table 5). 

Therefore, if we accept that such mass concentrates 
over a relaxation time, we can compare such estimates 
with the expected mass of a NSC as obtained by extrap¬ 


olation of the available NSC-galaxy host observed cor¬ 
relations. Indeed, many authors have provided scaling 
laws connecting the NSC mass with the host quantities 
(total mass, bulge mass, velocity dispersion) using high- 
resolution observations (Bdker et al. 2002; Cote et al. 
2006; Ferrarese et al. 2006; Bekki & Graham 2010; Seth 
et al. 2010; Leigh et al. 2012; Turner et al. 2012; Erwin & 
Gadotti 2012; Scott & Graham 2013), or theoretical ar¬ 
guments (Leigh et al. 2012; Antonini 2013; Arca-Sedda & 
Capuzzo-Dolcetta 2014b). Using the Leigh et al. (2012) 
correlation between the NSC mass and the host galaxy 
mass: 

Log(MNsc/MQ) = 1.18Log(Mg/MQ) - 4.376, (11) 

and the theoretical correlation found by Arca-Sedda & 
Capuzzo-Dolcetta (2014b): 

Log(MNsc/MQ) = 1.049Log(Mg/MQ) - 3.23, (12) 

and assuming for Henize 2-10 Mg = 6.4 x 10® Mq, we 
found Mnsc 10^ Mq using both the equations above. 
Under the hypothesis that the mass enclosed within 20 
pc can concentrate over a relaxation time, and adding 
the mass of the background galaxy contained within 10 
pc to the available mass in form of star cluster, we found 
a central mass of Mnsc ~ 8 x 10® Mq, a value in good 
agreement with that expected from scaling relations. 

5. GENERAL REMARKS AND GONGLUSIONS 

The aim of this work was the investigation of the pos¬ 
sible, future formation of a bright nucleus in the galaxy 
Henize 2-10 via orbital decay and merger of some of 
its stellar clusters. As an initial step, we used semi- 
analytical estimates, showing that the whole SCS of this 
galaxy should “collapse” to the galactic center in a time 
ranging between ^ 0.1 Gyr to ^ 1 Gyr, slightly depen¬ 
dent on the initial conditions chosen for the set of clus¬ 
ters. The predicted mass contribution of orbitally segre¬ 
gated clusters to the galactic nucleus has been estimated 
in the range 4 x 10® Mq to 6 x 10® Mq dependending on 
the initial conditions. 

To test quantitativaley these prediction we performed 
direct summation A-body simulations. These simula¬ 
tions constitute a direct test of the dry-merger scenario, 
and show that the decay and merging of star clusters 
would actually give rise to a detectable, bright nucleus 
whose mass is compatible with the observed masses of 
NSGs. 

The aim of the paper was also to understand whether 
and how the presence of a central massive BH would 
affect the formation process of a NSG in Henize 2-10. 
The role of a BH in the galactic center has been 
already studied and shown to have a relevance on the 
infalling clusters (see for example Antonini et al. (2012); 
Antonini & Merritt (2012); Arca-Sedda & Capuzzo- 
Dolcetta (2014a); Antonini (2014)). The effect consists is 
both a modulation of the dynamical friction on orbiting 
cluters and an increased efficiency of the tidal erosion 
process, leading (in some cases) to disruption of clusters 
before they reach the galactic center (Tremaine et al. 
1975; Antonini 2013; Arca-Sedda & Capuzzo-Dolcetta 
2014b) thus reducing or even preventing the formation 
of a NSC therein. 

In one of the simulations performed here (SI), we 
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Figure 14. As in Figure 3, but in the case of configuration S3. In this model the set of initial conditions for the SSCs is the same as in 
configuration SI, but the galaxy model does not host a MBH at its center. As consequence, the clusters reach the galactic center with a 
fraction of their initial mass larger than the estimates obtained from configuration SI. 


showed that it is possible to form a quite evident, bright, 
nucleus in the Henize 2-10 center within about 100 Myr. 
Since observations seem to suggest that the Henize 2-10 
clusters he on the same plane, we investigated such 
a possibility in simulation S2. In this case, the mass 
deposition to the central galactic region is significantly 
larger than in configuration SI, due to the smaller 
apocenters of the SSCs orbits. 

Once the NSC has been formed around the galac¬ 
tic center, its secular evolution is mainly driven by 
two-body relaxation, which brings to a further central 
concentration of the nucleus. This might lead to the 
formation of a very dense nucleus and, perhaps, to a 


slow accretion of stars onto the MBH. 

Our main results can be summarized as follows: 

• We created a model of the inner region of the 
galaxy Henize 2-10 using the highest resolution 
data available, representing this central region with 
more than 10® particles, which guarantees a high 
level in our A^-body modelization. 

• Our direct summation Al-body simulations show 
that the future evolution of the SCS in such galaxy 
will likely lead to the formation of a bright nucleus 
on a time-scale smaller than 1 Gyr. The growth of 
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Figure 15. Cumulative projected mass distribution of the SCS 
at different times in the configuration S3. The x-axis reports the 
projected distance from the galactic center. Most of the deposited 
mass is packed within i? = 20 pc from the cluster center within 
~ 40 Myr, while at later time such mass concentrates toward the 
innermost region {R < 10 pc). 



Figure 16. The same as in Figure 8, but for the configuration S3. 



Figure 17. Initial (straight line) and final (dashed line) surface 
density profile of the system galaxy+SCS in the simulation S3. 
The x-axis is the projected distance from the galactic center. It 
is clearly visible an overdensity that extends up to ~ 6 pc, with a 
total mass of M^gc ~ 5.1 X 10° Mq and half-light radius rpjgc ~ 2 
pc. 


the nucleus is maximized if the clusters lie, initially, 
on the same plane with nearly circular orbits. 

• The major contribution to the NSC formation is 
given by the merger of heaviest clusters (C1-C4), 
while lighter clusters experience a strong mass loss 
that leads to an almost complete disruption before 
they can reach the BH. 

• The decay and merging of massive clusters leads 
to a clearly detectable nucleus both in simulations 
SI and S2, whose projected density profiles extends 
out to 10 pc from the BH and whose effective radius 
is in the range 2 — 5 pc. Intriguingly, we found that 
the bright nucleus has a mass comparable to the 
BH mass. 

• The “future” NSC in simulations SI and S2 will 
have a mass in agreement with the estimates given 
by observational and (approximated) theoretical 
estimates. 

• Motivated by the uncertainty in the evaluation of 
the central BH mass, we discussed the role of a 
variation of the mass of the BH candidate in Henize 
2-10 galaxy, highlighting that, in the observed lim¬ 
its on the BH mass, we do not expect a significant 
change in the formation process of a NSC. In con¬ 
figuration S3, we used the same set of initial condi¬ 
tions as in configuration SI assuming, in this case, 
that Henize 2-10 does not host a MBH at its cen¬ 
ter. We found that the presence of a MBH slightly 
affects the decay process, since the full orbital de¬ 
cay occurs over a time which is almost the same 
in the two configurations, at least for the most 
massive clusters. On the other hand, our results 
showed that the MBH contributes significantly to 
the tidal erosion of the heaviest clusters that, in 
configuration SI, reach the galactic center with a 
mass ~ 15% smaller than in configuration S3. 

• Finally, we found that the mass-loss process for the 
lighter clusters is mostly determined by the tidal 
forces of the galactic background. Indeed, lighter 
clusters are efficiently disrupted even in the con¬ 
figuration S3. This is likely due to the fact that 
clusters C5-C11 spend most of their time on or¬ 
bits that never get very close to the galactic cen¬ 
ter, where the tidal effect arising from the MBH 
is maximized, and therefore the disruption process 
is mostly due to the tidal forces induced by the 
galactic background. 

6 . ACKNOWLEDGEMENTS 

MAS acknowledges the support provided by MIUR 
through the grant PRIN 2010 LY5N2T 005. 

REFERENCES 

Aharon D., Perets H. B., 2015, ApJ, 799, 185 1 
Antonini F., 2013, ApJ, 763, 62 1, 3.4, 4, 5 
Antonini F., 2014, ApJ, 794, 106 5 

Antonini F., Capuzzo-Dolcetta R., Mastrobuono-Battisti A., 
Merritt D., 2012, ApJ, 750, 111 1, 5 
Antonini F., Merritt D., 2012, ApJ, 745, 83 5 
Arca-Sedda M., Capuzzo-Dolcetta R., 2014a, ApJ, 785, 51 2.3, 
2.3, 3.2, 4, 3.4, 5 


































14 


Arca-Sedda, Capuzzo-Dolcetta, Antonini, Seth 


Arca-Sedda M., Capuzzo-Dolcetta R., 2014b, MNRAS, 444, 3738 
1, 3.4, 4, 4, 5 

Bekki K., Couch W. J., Shioya Y., 2006, ApJ, 642, L133 1 
Bekki K., Graham A. W., 2010, ApJ, 714, L313 4 
Binney J., Tremaine S., 2008, Galactic Dynamics: Second 
Edition. Princeton University Press 3.1 
Boker T., Laine S., van der Marel R. P., Sarzi M., Rix H.-W., Ho 
L. C., Shields J. C., 2002, AJ, 123, 1389 1, 4 
Buonanno R., Corsi C. E., Castellani M., Marconi G., Fusi Pecci 
F., Zinn R., 1999, AJ, 118, 1671 1 
Buonanno R., Corsi C. E., Zinn R., Pecci F. F., Hardy E., 
Suntzeff N. B., 1998, ApJ, 501, L33 1 
Capuzzo-Dolcetta R., 1993, ApJ, 415, 616 1, 2.3, 3.4 
Capuzzo-Dolcetta R., Mastrobuono-Battisti A., 2009, A&A, 507, 
183 1 

Capuzzo-Dolcetta R., Miocchi P., 2008a, ApJ, 681, 1136 1 
Capuzzo-Dolcetta R., Miocchi P., 2008b, MNRAS, 388, L69 1 
Capuzzo-Dolcetta R., Spera M., Punzo D., 2013, Journal of 
Computational Physics, 236, 580 3.1 
Capuzzo-Dolcetta R., Vicari A., 2005, MNRAS, 356, 899 2.3 
Chandar R., Leitherer C., Tremonti C., Calzetti D., 2003, ApJ, 
586, 939 2.3 

Chandrasekhar S., 1943, ApJ, 97, 255 2.3 
Cole D. R., Dehnen W., Read J. L, Wilkinson M. I., 2012, 
MNRAS, 426, 601 1 

Coleman M. G., Da Costa G. S., Bland-Hawthorn J., Freeman 
K. G., 2005, AJ, 129, 1443 1 

Cote P., Piatek S., Ferrarese L., Jordan A., Merritt D., Peng 
E. W., Ha§egan M., Blakeslee J. P., Mei S., West M. J., 
Milosavljevic M., Tonry J. L., 2006, ApJS, 165, 57 1, 4 
Cresci G., Vanzi L., Sauvage M., Santangelo G., van der Werf P., 
2010, A,^A, 520, A82 3.3 
Dehnen W., 1993, MNRAS, 265, 250 2.1 

den Brok M., Peletier R. F., Seth A., Balcells M., Dominguez L., 
Graham A. W., Carter D., Erwin P., Ferguson H. C., et al . 
2014, MNRAS, 445, 2385 1 

Dinescu D. I., Keeney B. A., Majewski S. R., Girard T. M., 2004, 
AJ, 128, 687 1 

Erwin P., Gadotti D. A., 2012, Advances in Astronomy, 2012 1, 4 
Ferrarese L., Cote P., Dalla Bonta E., Peng E. W., Merritt D., 
Jordan A., Blakeslee J. P., Ha§egan M., Mei S., Piatek S., 
Tonry J. L., West M. J., 2006, ApJ, 644, L21 1, 4 
Goerdt T., Moore B., Read J. I., Stadel J., Zemp M., 2006, 
MNRAS, 368, 1073 1 

Graham A. W., Spitler L. R., 2009, MNRAS, 397, 2148 1 
Hopkins P. F., Quataert E., 2010, MNRAS, 407, 1529 1 
Johnson K. E., Kobulnicky H. A., 2003, ApJ, 597, 923 1 
Johnson K. E., Leitherer C., Vacca W. D., Conti P. S., 2000, AJ, 
120, 1273 1 

King A., 2003, ApJ, 596, L27 1 
King A., 2005, ApJ, 635, L121 1 
King I. R., 1966, AJ, 71, 64 2.2 

Kobulnicky H. A., Dickey J. M., Sargent A. I., Hogg D. E., Conti 
P. S., 1995, AJ, 110, 116 1, 3.3 
Kobulnicky H. A., Martin C. L., 2010, ApJ, 718, 724 1 


Kormendy J., Ho L. C., 2013, ARA&;A, 51, 511 1 
Leigh N., Boker T., Knigge C., 2012, MNRAS, 424, 2130 1, 4 
Mackey A. D., Gilmore G. F., 2003, MNRAS, 340, 175 1 
Marquart T., Fathi K., Ostlin G., Bergvall N., Gumming R. J., 
Amram P., 2007, A&A, 474, L9 2.1 
Merloni A., Heinz S., di Matteo T., 2003, MNRAS, 345, 1057 1 
Milosavljevic M., 2004, ApJ, 605, L13 1 

Miocchi P., Lanzoni B., Ferraro F. R., Dalessandro E., Vesperini 
E., Pasquato M., Beccari G., Pallanca C., Sanna N., 2013, ApJ, 
774, 151 2.2 

Nayakshin S., Wilkinson M. L, King A., 2009, MNRAS, 398, L54 
1 

Neumayer N., Walcher C. J., 2012, Advances in Astronomy, 2012 
1 

Nguyen D. D., Seth A. C., Reines A. E., den Brok M., Sand D., 
McLeod B., 2014, ApJ, 794, 34 1, 2.1, 2.1, 2.2, 2.2, 3.3, 3.4 
Ostriker J. P., Binney J., Saha P., 1989, MNRAS, 241, 849 2.3 
Ott J., Walter F., Brinks E., 2005a, MNRAS, 358, 1423 1 
Ott J., Walter F., Brinks E., 2005b, MNRAS, 358, 1453 1 
Perets H. B., Mastrobuono-Battisti A., 2014, ApJ, 784, L44 1 
Pesce E., Capuzzo-Dolcetta R., Vietri M., 1992, MNRAS, 254, 

466 2.3 

Portegies Zwart S. F., Baumgardt H., Hut P., Makino J., 
McMillan S. L. W., 2004, Nature, 428, 724 1 
Read J. L, Goerdt T., Moore B., Pontzen A. P., Stadel J., Lake 
G., 2006, MNRAS, 373, 1451 1 
Reines A. E., Deller A. T., 2012, ApJ, 750, L24 1 
Reines A. E., Sivakoff G. R., Johnson K. E., Brogan C. L., 2011. 
Nature, 470, 66 1, 3.4 

Richstone D., Ajhar E. A., Bender R., Bower G., Dressier A., 
Faber S. M., Filippenko A. V., Gebhardt K., Green R., Ho 
L. C., Kormendy J., Lauer T. R., Magorrian J., Tremaine S., 
1998, Nature, 395, A14 1 

Rossa J., van der Marel R. P., Boker T., Gerssen J., Ho L. C., 

Rix H.-W., Shields J. C., Walcher C.-J., 2006, AJ, 132, 1074 1 
Santangelo G., Testi L., Gregorini L., Leurini S., Vanzi L., 
Walmsley C. M., Wilner D. J., 2009, A&A, 501, 495 1, 2.2 
Scott N., Graham A. W., 2013, ApJ, 763, 76 1, 4 
Seth A., Agiieros M., Lee D., Basu-Zych A., 2008, ApJ, 678, 116 
1, 2.3 

Seth A., Cappellari M., Neumayer N., Caldwell N., Bastian N., 
Olsen K., Blum R., Debattista V. P., McDermid R., Puzia T., 
Stephens A., 2010, in Debattista V. P., Popescu C. C., eds, 
American Institute of Physics Conference Series Vol. 1240 of 
American Institute of Physics Conference Series, Nuclear Star 
Clusters and Black Holes, pp 227-230 4 
Spera M., Arca-Sedda M., Capuzzo-Dolcetta R., 2015, ArXiv 
e-prints 3.4 

Tremaine S. D., 1976, ApJ, 203, 72 1 

Tremaine S. D., Ostriker J. P., Spitzer Jr. L., 1975, ApJ, 196, 407 
1, 2.3, 5 

Turner M. L., Cote P., Ferrarese L., Jordan A., Blakeslee J. P., 
Mei S., Peng E. W., West M. J., 2012, ApJS, 203, 5 1, 4 
Vacca W. D., Conti P. S., 1992, ApJ, 401, 543 2.1, 2.1 


